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We study the phenomenon of duality in hard exclusive reactions to which QCD factorization 
^ ■ applies. Considering "two-photon" -like processes in the scalar ip% model and also two-hadron (pion) 

f*^ ' production from the collisions of a real (transversely polarized) and a highly virtual, longitudinally 

i polarized, photon in QCD, we identify two regimes of factorization each of them associated with 

a distinct nonperturbative mechanism. One mechanism involves twist-3 Generalized Distribution 
Amplitudes, whereas the other one employs leading-twist Transition Distribution Amplitudes. In 
the case of the scalar ip% model, we find duality in that kinematical region where the two mechanisms 
overlap. In the QCD case, the appearance of duality is sensitive to the particular nonperturbative 
■ model applied and can, therefore, be used as an additional adjudicator. 

PACS numbers: 13.40.-f,12.38.Bx,12.38.Lg 
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I. INTRODUCTION 
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in: 

Lacking a complete theoretical understanding of color confinement, the only method of applying QCD is based 
| on the factorization of the dynamics and the isolation of a short-distance part which is amenable to perturbative 

■ techniques of quantum field theory (see, [l|, 5 S] and for a review, for instance, Q and references cited therein). 
On the other hand, the long-distance part has to be parametrized in terms of matrix elements of quark and gluon 

■ operators between hadronic states (or the vacuum). These matrix elements have nonperturbative origin and 
have to be either extracted from experiment or be determined by lattice simulations. In many phenomenological 
applications they are usually modeled by applying various nonperturbative methods. 

The QCD description of more involved hadronic processes requires the introduction of new hard parton sub- 
processes and nonperturbative functions. Important particular examples are exclusive hadronic processes which 
involve hadron distribution amplitudes (DAs), generalized distribution amplitudes (GDAs), and generalized 
parton distributions (GPDs) [ElElzL IEIEI- Within this context, collisions of a real and a highly- virtual photon 
provide a useful tool for studying a variety of fundamental aspects of QCD. A special case of this large class of 
processes is the exclusive two-hadron (pion) production in that region, where one initial photon is far off-shell 
(with its virtuality being denoted by Q 2 ) and longitudinally polarized, with the other photon being transversely 
polarized, while the overall energy (or, equivalently, the invariant mass of the two hadrons), s, is small. Such 
a process factorizes into a perturbatively calculable short-distance dominated subprocess, that describes the 
scattering 7*7 — * qq or 7*7 — > gg, and nonperturbative matrix elements measuring the transitions qq — > AB 
and gg — > AB. 

Another extensively studied — both theoretically and experimentally — two-photon process is the pion-photon 
transition form factor. This exclusive process has been measured by the CELLO [l(| and the CLEO [Tlj 
Collaborations for an almost on-shell (i.e., real) photon, whereas the other one has a large virtuality of up 
to 9 GeV 2 and is transversely polarized. These data have been analyzed over the years by many authors 
using various theoretical approaches — see, for instance, [U, EE Q El EE E3, El, EE HE HH- The major 
outcome of the most recent study [2(| is that the CLEO data exclude the Chernyak-Zhitnitsky pion distribution 



amplitude 22| at the 4<r level, while, perhaps somewhat surprisingly, also the asymptotic distribution amplitude 
seems to be disfavored (being off at the 3tr level). These findings are supported by the latest high-precision 
lattice simulations by two independent collaborations [IE [24[ • It is remarkable that the only pion distribution 
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FIG. 1: (Top) Dual factorization mechanisms in describing the fusion of a real photon with a highly-virtual and longi- 
tudinally polarized photon in the Euclidean ip 3 E model. (Bottom) The analogous situation in QCD. 



amplitude found to be within the la error ellipse of the CLEO data, while complying: with the new lattice 
results just mentioned, is the pion model derived from nonlocal QCD sum rules [25|. 

These results lend credibility to the usefulness of two-photon processes as a means of accessing and under- 
standing the parton structure of hadrons. On that basis, one may try to assess more complicated two-photon 
processes with two pions in the final state. Indeed, recently, nonperturbative quantities of a new kind were 
introduced — transition distribution amplitudes (TDAs) [2^, [23, [2^] — which are closely related to the GPDs and 
describe the transition (77 — ► qA. In contrast to the GDAs, mentioned above, such nonperturbative objects, 
like the TDAs, appear in the factorization procedure when the Mandelstam variable s is of the same order of 
magnitude as the large photon virtuality Q 2 , while t is rather small. Actually, there exists a reaction where both 
amplitude types, GDAs and TDAs, can overlap. Precisely this can happen in the fusion of a real transversely 
polarized photon with a highly- virtual longitudinally polarized photon that gives rise to a final state comprising 
a pair of pions. This reaction can potentially follow either path, i.e., it can proceed via the twist-3 GDAs, or 
through the leading-twist TDAs, as illustrated in Fig. [T] The first mentioned possibility is linked to the case 
when one of the involved photons is longitudinally polarized, whereas the other one is transversely polarized. 
The second possibility corresponds to the same situation, but appears in contrast to that in leading- twist order, 
as we will show in Sec. IIIII Thus, a comparative analysis of both possible mechanisms for this reaction seems 
extremely interesting both theoretically and phenomenologically. 

A simultaneous analysis of these two mechanisms in the production of a vector-meson pair was carried out 
in [29| . The authors found that these mechanisms can be selected by means of the different polarizations of the 
initial-state photon. In the case of (pseudo)scalar particles, like the pions, this effect is absent and, therefore, 
this opens up a window of opportunity to access the overlap region of both mechanisms and their duality as 
opposed to their additivity. By duality we mean here that adding the contributions of the two mechanisms 
would lead to a double-counting. Hence, these mechanisms are in antagonism to each other and should be 
considered as two alternative ways to describe the same physics. 

In this paper, we will verify the possibility for duality between the two different mechanisms of factorization, 
associated either with GDAs, or with TDAs, in the regime where both Mandelstam variables s and t are rather 
small compared to the large photon virtuality Q 2 . This will be done below by considering first (Sec. [TTJ) the 
Euclidean (^-analogue of QCD, followed then by an investigation of the exclusive pion-pair production in 77^ 
collisions in Sec. IIIII Our findings are summarized and further discussed in Sec. IIVI where we also provide our 
conclusions. 
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II. REGIMES OF FACTORIZATION WITHIN THE v^-MODEL 

We start our analysis with the factorization of a two-photon process in which the two photons undergo a 
fusion with the subsequent creation of two pions. For the sake of simplicity, we work first in the scalar Euclidean 
ip E model, regarding it as a sort of a toy model in order to prepare the ground for the real QCD application to 
follow in the subsequent section. 

The most convenient way to analyze this four-particle amplitude in the ip E model is to work within the a- 
representation, following the lines of thought given in Q in the discussion of Deeply Virtual Compton Scattering. 
The reason is that within this representation (i) the calculations can be performed in a systematic way, not 
only at one loop but also at higher-loop levels, (ii) the factorization of the process can be studied in detail, and 
(iii) one can study the spectral properties of the nonperturbative input (GPDs, GDAs) because these properties 
are insensitive to the numerators of the quark and gluon propagators. In addition, complications owing to the 
spin structure do not affect them. It is worth recalling in this context that the a-representation of Feynman 
diagrams is not only a schematic means for identifying relevant integration regions in a factorization procedure, 
but can, in fact, serve to provide rigorous proofs of factorization theorems [![. 

The contribution of the leading "box" diagram, see Fig. [3J can be written as 
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where m 2 serves as a infrared (IR) regulator, s > 0, t > are the Mandelstam variables in the Euclidean region, 
and D = a\ + ... + a.4,. 





(a) (b) (c) 

FIG. 2: (a) The box diagram in a-space within the (^-theory. The left open lines in each graph simulate the "photons" 
with the appropriate kinematics (see text), whereas the open lines at the right denote the produced mock "pions". (b) 
"GDA factorization" indicated by an oval, (c) "TDA factorization" indicated by an oval. 



It is useful to introduce the dimensionless parameters 
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(here, the dimensional parameter A 2 is an arbitrary one) and, then, recast the a-parameters into the form 
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As a result, we get the following expression for the considered amplitude: 
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FIG. 3: The integrand of the exact amplitude (j4]) as a function of ft and ft in regime (a). Left panel: top view; right 
panel: 3D illustration. 



where D = ftftftft + ... + 0x020304 and ft = 1 — ft. The particular feature of the considered process is 
that the transferred momentum q 2 = Q 2 is large compared to the mass scale m 2 , the latter simulating here the 
typical scale of soft interactions. With respect to the other two kinematical variables s and t, one can identify 
four distinct regimes: 

(a) s <C Q 2 while t is of order Q 2 ; 

(b) t <^Q 2 while s is of order Q 2 ; 

(c) s, t < Q 2 ; 

(d) s, t are both of order Q 2 (this regime corresponds to large-angle scattering Q and will not be addressed 
here) . 

Notice that within the regimes (a) and (b), amplitude JT]) or Q can be factorized. Recall that the factorization 
of ((T|) and (|3]) at the leading-twist level in the a- or /^-representation is equivalent to the calculation of the 
leading Q 2 -asymptotics of the Laplace-type integral 



oo 

F(X) = J da ff («)exp [-Xf(a)} 
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with a large and positive parameter A and the function /(a) having a minimum at the point a — 0. 
Regime (a) 

Let us begin the discussion of this regime by considering how the integrand in ([4]) depends on the parameters 
ft and ft, keeping in mind that the parameters ft and ft will be integrated out. To model this regime, we fix, 
for example, the value of s/A 2 to be equal to 1 and set t/A 2 = 50.0, Q 2 /A 2 — 100.0. Then, the two-variable 
function representing the integrand in (|4|) behaves as shown in Fig. [2l From this figure, one can conclude that 
the main contribution to the exact amplitude (j4|) emanates from the area restricted by small ft values around 
Pi ~ 0.05 and ft values varying in the wide interval [0, 0.9]. We will call this support of the integrand ft-wing. 
Indeed, we may divide the whole support into the following domains (see Fig. [4]) 

Region- 1: ft e [0, 0.1], ft G [0, 0.1]; 
Region - 2 : ft e [0.1, 1], ft G [0.1, 1]; 
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Region - 3 : (3 X € [0, 0.1], #> € [0.1, 1]; 
Region - 4 : ft£ [0.1, 1], (3 2 € [0, 0.1] 



(5) 



Our numerical computation shows that the integrations over the strip including the regions 1 and 3 (the /So- 
wing) gives a 85% contribution to the amplitude ([U relative to the contribution from the whole support. Thus, 
we may approximate the exact amplitude by integrating over the /?2-wing only. [Note that the contribution 
from the integration over the region 1 provides approximately a mere 5 %.] 

ft 



0.1 



0.1 



1 p, 



FIG. 4: Regions of integration in f3 space of the integrand of the exact amplitude (|4| as a function of /3i and ,02. 

On the other hand, in this regime, the exact amplitude can also be estimated by the asymptotic formula which 
results from the analytic integration over the region where (3\ ~ (corresponding to the "GDA factorization" — 
recall Fig. [J): 
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where D = (3 2 (3~3(3~4 + (3 2 f3 3 (3 4 + (3 2 (3s,f34 L . Indeed, the integration over (3\ ~ eliminates the Q 2 -dependence in 
the exponential function and provides the main contribution to the amplitude A. Schematically, this means 
that the propagator parameterized by ot\ or (3\ can be associated with the partonic (hard) subprocesses, while 
the remaining propagator constitutes the soft part of the considered amplitude, i.e., the scalar version of the 
GDA (see Fig.[2jb)). Note that the considered process is going through the s-channel, which would correspond 
to the GDA mechanism in QCD. 1 

Concluding this discussion, one should stress that when s/A 2 = 1.0, i/A 2 = 50.0, the ratio between the 
asymptotic ([6]) and the exact amplitude (|4]) is 



R = 



A 



GgA =1.01. 

A 



(7) 



This ratio shows that the asymptotic formula ^ reproduces the exact amplitude ^ with a rather high accuracy, 
though slightly overestimating it. Such a situation is very natural for the considered regime of the amplitude 
within the a parametrization and is not at odds with other factorization procedures. 
Regime (b) 

Here, we have to eliminate from the exponential in Eq. ^ the variables Q 2 and s, which are large. This can 



In anticipation of the QCD case later on, we already use the notions GDA and TDA, though, strictly speaking, they are not 
applicable here. 
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FIG. 5: The integrand of the exact amplitude (j4| as a function of 0\ and /3;2 in regime (c). 



be achieved by integrating over the region 02 ~ 0, i.e., by cutting the line corresponding to this parameter. 
Performing similar manipulations as in regime (a), we find that the scalar TDA amplitude (associated with the 
factorization in the t-channel (see Fig. Hfc))) can be related to the scalar GDA via 

i? DA (s, t, m 2 ) = A^ DA (t : s, m 2 ) . (8) 

Regime (c) 

The most important regime for our investigations on duality is when it happens that both variables s and t are 
simultaneously small compared to Q 2 (but still much larger than the soft scale m 2 ): s, t -C Q 2 . The question 
which may arise in this case is what happens with factorization. In other words, can the asymptotic formula be 
useful even in this case? 

In investigating this issue, we face two different options. 

Option 1 

Start with, say, the regime (a) , since there is no doubt that the asymptotic formula (or the factorization) can be 
applied to both regimes (a) and (b) . We begin by considering the integrand of @ in terms of the two variables 
0i and 02- Then, we tune the value of the variable t/A 2 down to small enough values in order to understand the 
behavior of the integrand in this regime. We originally designed the integrand of for t/A 2 assuming values 
within the set {25, 10, 5, 2, 1}, but for our purposes it is necessary to study the special case of t/A 2 = s/A 2 = 1 
presented in Fig. El From this graphics we see that the support of the considered function contains now two 
wings from which the function becomes saturated. Note that the 0\- and /^-wings are the same in this case 
and that both wings are equally essential. The numerical analysis shows that the integration in (|4]) over the 
regions 1, 3, and 4 (0i and 02 wings), i.e., 

Reg.-l Rcg.-3 Reg. -4 

[ Q 2 PlfoMi S 02040103 t 01030204, m 2 D 

CXP [ A 2 D A 2 D A 2 f) A 2 fcfofafa. 

provides approximately a 80% contribution to the exact amplitude in this regime. Therefore, the saturation of 
the exact amplitude ^ in regime (c) results from contributions originating from both wings 0\ and 02- Hence, 
restricting attention to only one of them is not permissible. In this sense, we deal with a particular case of 
"additivity" of the two wings. 
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Option 2 

However, there exists an alternative way to determine the exact amplitude, based on the analytic integration 
of Q in /3 space. Namely, we may continue the asymptotic formula ([5]), obtained within regime (a), to regime 
(c), where t/A 2 and s/A 2 are the same (or approximately the same). Then, we will observe that the asymptotic 
formula ^ still describes the amplitude 



A « A^ DA (s,t,m 2 



(10) 



with an acceptable accuracy, still slightly overestimating the exact amplitude ((3]). This becomes evident by 
analyzing the ratio ([7]) in terms of t/A 2 , as depicted in Fig. [5J One sees that the large i-tail of ratio ([7]) is 
actually well-described by the findings in regime (a). At the same time, when we go to the regime where t is 
small or of the same order as s, the asymptotic formula ([6j> — which was continued from regime (a) — describes 
even in this case the exact amplitude with a tolerable 10% deviation. Indeed, in the small t region, the ratio 
([7]) increases from approximately 1.01 to about 1.09 at t/A 2 w 1. Hence, we may conclude that, despite the 
fact that the exact amplitude ^ in regime (c) is saturated from contributions stemming from both wings /3i 
and /?2, the asymptotic formula ([6]), derived in regime (a) and continued to regime (c), is still suitable for an 
acceptable approximation of the exact amplitude. 

Due to the symmetry of regimes (a) and (b) under the exchange of the variables s <-> t, a similar statement is 
valid as regards the asymptotic formula derived in regime (b) — cf. ([8]). Therefore, in regime (c), the asymptotic 
expression ^ still remains a good approximation of the exact amplitude, so that one does not have to reestimate 
it using as an alternative Eq. ([5]). Insisting to do so, one would face a "double-counting" problem. Thus, in this 
sense, we deal with duality between © and © . Note that this duality is not a result of the symmetry between 
regimes (a) and (b). This symmetry merely simplifies the analysis of regime (b). 

Since the derivation of the asymptotic formula resembles the factorization procedure for the given amplitude 
(in which the asymptotic formula plays the role of the partonic part), the GDA and TDA factorizations within 
the (/3'g-model in regime (c) are equivalent to each other without a kinematical or dynamical prevalence of one 
over the other. This can be interpreted as a sort of duality between the GDA and TDA factorizations. 

In concluding this discussion, let us stress that there is a difference between the two possible ways of factoriza- 
tion, which are usually considered to be equivalent, because, there is either dominance of some parts of the full 
integration region, see @, or the application of the asymptotic (partonic) formula, see (fTU)) can (approximately) 
describe the exact process amplitude. In the case of regime (c), the first way of factorization implies additivity 
of the two wings. Alternatively, the continuation of the asymptotic (partonic) expressions gives rise to duality 
between them. 



III. TDA-FACTORIZATION VS. GDA- FACTORIZATION FOR 77* -> tvtv 

In the preceding section, we demonstrated how the duality between the GDA and the TDA-factorization 
works within a toy model. We now study the duality phenomenon in the case of real QCD. To this end, we 
consider the exclusive ir + ir~ production in a 7t7£ collision, where the virtual photon with a large virtuality Q 2 
is longitudinally polarized, whereas the other one is quasi real and transversely polarized. As already mentioned, 
the key feature of this process is that it exhibits two different kinds of factorization, based, respectively, on the 
GDA or the TDA mechanism (see Fig. [7]), with a potential overlap of both mechanisms in the kinematical 
regime where t and s are small. Notice that this situation differs from the cases studied in [5^1, where the 
additivity of the GDA and the TDA factorizations was explored and found that each factorization mechanism 
can be selected on account of the helicity of the particles. 

In contrast, in our case, the GDA and the TDA regimes correspond to the same helicity amplitudes, resembling 
the situation faced in the scalar case considered in the preceding section where there is only one helicity 
amplitude. To continue, let us recall the parametrization of nonperturbative matrix elements entering our 
analysis. Because the considered process involves a longitudinally and a transversally polarized photon, we 
have to deal with twist-3 parametrizing functions. To this end, we adhere to (30j . where the genuine twist-3 
contributions in the 77* — * tttt process were studied in detail (see, also, [Uj]). As regards the twist-2 contribution, 
related to the meson DA, we use the standard parametrization of the 7r + -to-vacuum matrix element which 
involves a bilocal axial- vector quark operator [l| . 

On the other hand, the 7 — > tt~ matrix elements, entering the TDA-factorized amplitude, can be parameter- 
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FIG. 6: The ratio 7?, see ([7]), as a function oft/A 2 



ized in the form (cf. [27 



— T 

(7r-(p2)|V(-*/2)7a[-*/2; z/2]^(z/2)| 7 (g', e')> = je^pa^i^, £, t) , 

Jot 

^-(p 2 )|^(-z/2)7 Q75 [-z/2;z/2]^(z/2)|7(q', £ ')> = ^4 ' *TP a Ai(x,{, t) , 

Jot 

where the symbol = stands for the Fourier transformation (with the appropriate measure) and 

i—g/2 



[~z/2;z/2] = Pexp 



ig I dx^A^{x)t a 

2/2 



(11) 



(12) 



denotes a path-ordered gauge link to ensure gauge invariance. In the following, we choose to work in the axial 
gauge which entails that the gauge link is equal to unity. In the nonperturbative amplitudes Vy and A\ 




LO GDA fact. 
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FIG. 7: Two ways of factorization: (Left panel) via the GDA mechanism and via the TDA mechanism (right panel). 
Only a single typical diagram for each case is shown for the purpose of illustration. 
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denote, respectively, vector and axial- vector TDAs, whereas the relative momentum is defined as P = (p 2 +q')/2, 
and A = p2 — q' . Notice that TDAs in the ERBL region, i.e., — £ < x < £, include the so-called D-term, where 
the pole contribution, related to the spinless resonance (this being the pion in our case), is included l32l . |33| . 
Finally, to normalize the axial- vector TDA, Ax, we adopt the philosophy of [27], [U, HH, [36|, H3, [H, [3!| and 
express Ai in terms of the axial- vector form factor measured in the weak decay ir — > Ivi'y, i.e., 

l 

J dxAx{x, f, t) = 2f 7T F A (t)/m^, (13) 
-1 



where f n = 0.131 GeV, = 0.140 GeV, and F A (0) w 0.012 |40|,|4l|. 

The next objects of interest in our considerations are the helicity amplitudes that are obtained from the usual 
amplitudes after multiplying them by the photon polarization vectors: 

A i0 , ±) =e^T^\ ^ = (0,^,-^,0), e^=(M=,0,0,^). (14) 



Discarding the A 2 ^ corrections, the only remaining contribution stems from the axial- vector one {A\) given by 
Eq. (fTTj) . Thus, the helicity amplitude associated with the TDA mechanism reads 

TDA g 'b) ■ A r TDA TDA 2m c f f j 0ir(y) f j 4 / t ±\ ( e " e rf \ 

-1 

where 

and where we have employed the 1-loop a s ((5 2 ) in the MS-scheme with Aqcd = 0.312 GeV for Nf = 3 [42] ]. 
Note that there is only a mild dependence on Aqcd. Choosing, for instance, a somewhat smaller value, say, 
around 0.200 GeV, would shift the TDA curves shown in Fig. [5] slightly downwards. 

On the other hand, the helicity amplitude which includes the twist-3 GDA can be written as (see, for example, 

MM) 

^ = ^^^ DA , ^ A = 2^±^(et + e l) Jdy d ^(y,C,W^f- l f), (17) 





where 

2t 

2C-l = /3cos^ ro , cos^ ro = w2 + Q2 -1 (18) 

and the partial derivative is defined by d$ — d/d(2C, — 1). Here the amplitude $1(2/, C> W 2 ) encodes the matrix 
element of the twist-2 operator. In deriving (|17p . we have used the Wandzura-Wilczek approximation for the 
twist-3 contribution and took into account that in our system (p 2 —pi) T = 2(p 2 — q') T ■ Besides, we would like to 
stress that the amplitudes (fl"5|) and (fT7|) are both of order 1/Q 2 , which offers room for duality. As pointed out 
above, duality between expressions (|15[) and (|17p may occur in that regime where both Mandelstam variables s 
and t are much smaller in comparison to the large photon virtuality Q 2 . In terms of the skewness parameters 
£ and £, this would imply that £ ~ 1 and ( ~ 0. 

To estimate the relative weight of the amplitudes with TDA or GDA contributions, we have to model these 
non-perturbative objects. 

Let us start with the TDAs. As a first step, we assume a factorizing ansatz for the i-dependence of the TDAs 
and write 

A 1 (x,Z,t)=2 — FA(t)A 1 (x,Z), (19) 
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where the f-indcpcndcnt function A±(x,£) is normalized to unity: 

l 

dxA 1 {x,£) = l. (20) 



-l 



To satisfy this condition, we introduce a TDA defined by 

l 

A non-norm/ i \ f 

A ± (x, 1) = N [ ' M= dxAl on - notm {x, 1) . (21) 

-l 

We now focus on the discussion of the i-independent TDAs in (fT9)) . Since we are mainly interested in TDAs 
in the Efremov-Radyushkin-Brodsky-Lepage (ERBL) [HQ region £ = 1, it is instructive to choose the following 
parametrization 



^non-norm ^ i) = (i _ ^) ^ + ai Cf'(x) + a 2 C" 2 J ' ^ {x) + a 4 Cf'(x)J , (22) 

where ai, a 2l £14 are free adjustable parameters, encoding nonperturbative input, and the standard notations for 
Gegenbauer polynomials are used. One appreciates that the TDA in the form of Eq. (j2"2"|) amounts to summing 
a D-term, i.e., the term with the coefficient ai, and meson-DA-like contributions. Indeed, would we eliminate 
the term with ai, we would obtain the standard parametrization for a meson DA. On the other hand, keeping 
only the term with ai .would reproduce the parametrization for the D-term Therefore, for our analysis, we 
suppose that a\ = do which is equal to —0.5 in lattice simulations [43|. 2 With respect to the parameters a 2 
and 04, we allow them to vary in quite broad intervals, notably, 0,2 € [0.3, 0.6] and 0,4 € [0.4, 0.8], that would 
cover vector-meson DAs with very different profiles at a normalization scale /j, 2 ~ 1 GeV 2 (see, for example, 
PH). Note that the profiles obtained with ansatz (|22]) bear remarkable resemblance to the TDAs calculated 
within a non-local chiral quark-soliton model in the recent work of Ref. [34[ . 

As regards the function $i(z, £), which denotes the corresponding GDA in (jTTJ) , we take recourse to the 
following model 0] 

$!(z,C,T^ 2 ) =W f zz(2z-l) (B w {W 2 )e lSa ^ + B 12 (W 2 )e tS ^ w ^ P 2 (cos9^ , (23) 

where the phase shift of the tttt scattering is defined by S (W = 0.8) w f and S 2 (W = 0.8) w 0.03?r [H,^. 
In ([23]) . the function i?io corresponding to L = does not contribute to (fl~7|) and, therefore, can be discarded. 
On the other hand, the function Bi 2 can be estimated by the method described in j48|. Since Bi 2 corresponds 
to the two-pion state with L = 2 and our GDAs probe only the isoscalar channel, we derive for B\ 2 at W 2 in 
the region of the /2-meson [I G (J PC ) = + (2 ++ )] mass the following expression: 

2 10 9f 2irv ffrMj Tf 2 

B 12 (W ) = - {M ^_ w2)2+ 2 V 2M 2 j (24) 
where // 2 = 0.056 GeV, M/ 2 = 1.275 GeV, Tf 2 = 0.185 GeV, and the decay constant gf 27r7T is defined by 




' 24-7T 

9f^ = \iTW T ^^^ T(/ 2 -> tttt) = 0.85r /2 . (25) 



In addition, we also model the function B i2 with the help of the most simple ansatz (see, [47j ) 

B 12 (Q)=f^R n , (26) 



2 The case of a\ = —4/Nf, which holds within the chiral quark-soliton model |44|| , is also included in our analysis. We found that 
there is no much difference in comparison with the lattice model. 
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FIG. 8: Helicity amplitudes J rTDA [cf. Eq. JT5}] and T WK j cf Eq gjj ag functions of q2 j using 0l 
lattice simulations [43 ] . The value of s/Q 2 varies in the interval [0.06, 0.3]. 



-0.5 found in 



where denotes the fraction of momentum carried by the quarks and antiquarks in the pion. In our numerical 
analysis below, we use for this fraction R n = 0.5, which is taken from the leading-order Gliick-Reya-Schienbein 
approximation (see, for example, (47|). 

Before continuing with the numerical analysis, let us pause for a moment to discuss the validity and self- 
consistency of the approximations made. Start with the pion DA: For simplicity, the asymptotic pion DA in 
Eq. (fT5]) was used. Using other pion DAs [2(|, the contribution from the TDA factorization may increase from 
about 6% for the BMS model [25| to approximately 56% for the CZ model [23]. Second, continue with the 
GDA: At the level of accuracy relevant for our arguments, we used for the twist-2 GDA <I>i the asymptotic 
form, the justification being provided by quark-hadron duality arguments given in [49j ] . Third, we assumed the 
Wandzura-Wilczek approximation, which relies upon the assumption of the smallness of the genuine twist-3 
contributions. This assumption is supported by measurements of the gi structure function. The normalization 
of the GDA is determined by the momentum conservation which is rigorously valid. By performing a crossing 
from GPDs (where momentum conservation applies) to the GDA channel, we used the improved expression 
([2"4"| accounting for meson dynamics. Turn to the TDA: The first ingredient is the normalization (fT5|) , which 
is fixed by the data for the pion weak decay (see for details [27| )■ Finally, for the form of the TDA, we used a 
parametrization [cf . ([22]) ] , which can accommodate different available models [H], [H, [36|, [37| , within the range 
of accuracy of the order of 10%, the latter pertaining to our whole numerical analysis. 

We calculated both functions F TJ: > A and f GBA and show the obtained results in Fig. [5] for ai — —0.5 (which 
corresponds to the lattice result of [43|]) and employing for simplicity the asymptotic pion DA. The dashed line 
corresponds to the function JF TDA , where we have adjusted the free parameters to a 2 = 0.6 and 04 = 0.8. The 
results, obtained for rather small values of these parameters, are displayed by the broken lines in the same 
figures. The dotted line denotes the function jF TDA with a-i = 0.5 and 04 = 0.6, whereas the dash-dotted line 
employs a 2 — 0.3 and 04 = 0.4. For comparison, we also include the results for the .F GDA . In that latter case, 
the dense-dotted line corresponds to the GDA amplitude, where the expression for B12 has been estimated via 
Eq. ([2~4"|) , while the solid line represents the simplest ansatz (j2"6")l with — 0.5. The particular values of the 
coefficients 0,2,0-4, we used, are within the fiducial intervals specified above in conjunction with Eq. ([22]) . If, 
for example, we would have used instead values of these parameters outside their fiducial intervals, then the 
outcome would change by orders of magnitude and duality would be badly violated. 
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From this figure one may infer that (i) in the case when the parameter B\i (which parametrizes the GDA 
contribution) is estimated with the aid of the Breit-Wigner formula, see (|24|) , and, while s, t <g; Q 2 , there is 
duality between the GDA and the TDA factorization mechanisms. Hence, the model for B12 which takes into 
account the corresponding resonances can be selected by duality, (ii) In the case of the simplest model for 
-B12, evaluated for R n — 0.5 (see, (|2"r?)) ). the duality between the GDA and the TDA factorization mechanisms 
reaches maximally the level of 35%, provided the axial- vector TDAs are normalized in terms of the axial form 
factor Fa- On that basis, one may conclude that the simplest model for B12 is rather inconsistent with duality 
and cannot be considered as a realistic option. 

IV. SUMMARY AND CONCLUSIONS 

The process amplitude for the fusion of a highly virtual, but small s photon (longitudinally polarized), with a 
quasi real one (with a transversal polarization) can be factorized in two different ways. The study of this process 
within a toy ip E model has shown the existence of duality between the factorization along the s-channel and 
the t-channel — which serve to simulate, respectively, the GDA and the TDA factorization mechanisms in QCD. 
This observation persists also in the case of such reactions in real QCD. Indeed, the amplitude for the process 
77* — > itit can be factorized in two distinct ways, having recourse to the large virtuality of one of the photons 
involved, as compared to the scales of the soft interactions. As mentioned before, depending on the variables 
s and t, two different mechanisms for this reaction were identified. The first factorization mechanism, which 
employs GDAs, takes place when s<0 2 , while t is of the order of Q 2 . On the other side, when it happens that 
t <C Q 2 and s ~ Q 2 , the TDA mechanism of factorization becomes dominant in such a 77* collision. We have 
shown that when both variables s and t are simultaneously much smaller than Q 2 in the collision of a real and 
transversely polarized photon with a highly virtual longitudinally polarized photon in the pion-pair production, 
then both types of factorization are possible giving rise to duality in this region. 

Our observations in QCD can be summarized and outlined as follows. 

We have shown that when it happens that both Mandelstam variables s and t are much less than the large 
momentum scale Q 2 , with the variables s/Q 2 and t/Q 2 varying in the interval (0.001, 0.7), then the TDA and 
the GDA factorization mechanisms are equivalent to each other and operate in parallel. This marks a crucial 
difference between our approach and previous ones [2!| which dealt with the additivity of the TDA and the GDA 
factorization mechanisms. Note that the presented results were obtained using rigorous constraints and available 
nonperturbative models for the DAs, GDAs, and the TDAs. We found that duality appears only for values 
of the model parameters within their fiducial intervals, associated with each particular model. Therefore, we 
demonstrated that duality may serve as a tool for selecting suitable models for the non-perturbative ingredients 
of QCD factorization of various exclusive amplitudes. 

We also observed that twist-3 GDAs appear to be dual to the convolutions of leading-twist TDAs and DAs, 
multiplied by a QCD effective coupling — the latter stemming from the hard-gluon exchange. This reflects at 
the level of the particular process, we are considering, the manifestation of duality for the GPDs and GDAs 
themselves [H, Hfl |51| . One may ask to which extent the observed duality between different QCD factorization 
mechanisms is an exceptional finding or rather represents a generic feature of QCD dynamics. In this connection, 
one may recall the different ways of describing single-spin asymmetries (SSA) using either twist-3 collinear 
factorization 52, 53] or by employing the transverse-momentum dependent (TMD) Sivers distribution function 
(53 | . Because of the T-odd nature of the SSA and the Sivers function, the latter may emerge only due to the 
interactions between the hard and the soft par ts (55| of the process and represents an effective distribution [5(| . 
Moreover, its transverse moment is related [57] to a particular twist-3 matrix element. This sort of relation may 
be explored [58| as a matching procedure (or, in our language, as duality) between TMD factorization at low 
Pt-i on one hand, and collinear twist-3 factorization at large Pt, on the other. Recently [59], such a matching 
procedure was carried out in detail for various observables and it was found that duality strongly depends on 
the type of the considered observable. In the case of the Sivers function, there is a possibility of enlarging the 
duality intervals by either using the twist-3 approach at low Pt [13] or by applying the Sivers function at large 
Pt [6l|. The latter procedure finds justification in the framework of the twist-3 factorization, provided only the 
first non-vanishing transverse moment of the Sivers function is retained and that color factors, responsible for 
the effective nature of the Sivers function, are properly taken into account. 

Duality also emerges when another type of QCD factorization in semi-inclusive processes is considered, which 
makes use of so-called fracture functions. The momentum conservation guarantees that the description in 
terms of either fragmentation or fracture functions is complete and one should have duality between these two 
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mechanisms in the overlap region (see [62J and references therein). Bottom line: We have shown that duality 
leads to new insights into a number of phenomena in QCD, in particular to those related to exclusive processes 
and spin-dependent quantities. 
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